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Abstract: We examine some conditions which are needed to allow the Peccei-Quinn scale 
(f a ) as high as the SUSY GUT scale (M GUT ~ 10 16 GeV) in the context of the PQMSSM 
with either an axino or a neutralino LSP. The main problem in non-SUSY models with 
f a ~ Mqut is the generation of an overabundance of axion dark matter (~ fV & ) due to 
vacuum misalignment. We show that once all the components of the axion supermultiplet 
are included, the upper limit on f a can be evaded due to large entropy injection from 
saxion decays. This large entropy injection also dilutes all other quasi-stable relic densities, 
naturally evading the BBN constraints and solving the gravitino problem. We find that 
f a ~ Mqut can be allowed by relic density/BBN constraints provided that the saxion 
mass m s > 50 TeV, the initial saxion field value is of order of the PQ scale and the initial 
axion mis-alignment angle 0j < 0.05. These restrictions can be considerably loosened for 
f a ~ 10 15 GeV. The allowed range for the re-heat temperature (Tr) is strongly dependent 
on the nature of the LSP. For the axino LSP, Tr > 10 8 GeV, while for neutralino LSP any 
value is allowed. In the latter case, /„ ~ Mqut can be more easily accommodated than 
in the axino LSP scenario. For f a ~ Mqut in SUSY models, the dark matter abundance 
should be dominated by axions, albeit with mass m a ~ 10" 10 eV, far below the region 
currently probed by experiment. 
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1. Introduction 



The strong CP problem remains one of the central puzzles of QCD which evades ex- 
planation within the context of the Standard Model. The crux of the problem is that an 
additional CP violating term in the QCD Lagrangian of the form 1 6g 2 /32-K 2 G A "GA f iu ought 
to be present as a result of the t'Hooft resolution of the U{1)a problem via instantons and 
the vacuum of QCD [1]. Here, 6 = 9 + arg detAi actually consists of two terms: one 
from QCD and one from the electroweak quark mass matrix. The experimental limits on 
the neutron electric dipole moment however constrain \9\ < 10~ 10 [2]. Explaining why the 
sum of these two terms should be so small is the essence of the strong CP problem [3]. 

An extremely compelling solution proposed by Peccei and Quinn [4] is to hypothesize 
an additional global ?7(1)pq symmetry, which is broken at some mass scale f a > 10 9 GeV. 
A consequence of the broken PQ symmetry is the existence of a pseudo-Goldstone boson 
field: the axion a(x) [5]. The low energy Lagrangian then includes the interaction term 



a s a(x) 



£3 5Z77J7^^, (1.1) 



where a s = g 2 /Air, g s is the QCD coupling constant and N is the model-dependent color 
anomaly factor. (From here on we assume N = 1, but all our results can be extended 
to any ./V value with the replacement f a —> f a /N.) Since a(x) is dynamical, the entire 
CP-violating term settles to its minimum at zero, thus resolving the strong CP problem. 
A consequence of this very elegant mechanism is that a physical axion field should exist, 
with axion excitations of mass [6] 

ma = 6eV "^. (1 . 2) 

J a 

Due to its tiny mass, the interactions of the axion field need to be strongly suppressed- 
fa ^ 10 9 GeV- otherwise they would have a strong impact on low energy physics and 
astrophysics [7], such as too rapid cooling of stars and supernovae. On the other hand, 
under some model assumptions (see Sec. |3.1| below), the axion relic density Vt a h 2 < 0.1123 
bound requires f a < 10 12 GeV. As a consequence a new physics scale much larger than 
the electroweak scale (~ my/) nas to be introduced into the Standard Model [8,9]. This 
results in large radiative corrections to the Higgs mass, which then requires a large amount 
of fine-tuning to stabilize the EW scale. This leads to the well known hierarchy problem. 

So far, one of the most compelling solutions to naturally stabilize the EW scale is 
supersymmetry (SUSY), which reduces quadratic divergences to merely logarithmic, and 
ameliorates the fine-tuning problem [10]. The supersymmetric version of the Standard 
Model (MSSM) also has other compelling features, such as viable Dark Matter (DM) 
candidate(s) and unification of the gauge coupling constants at Mqut - 2 x 10 16 GeV. 
In order to accomodate the PQ solution in the MSSM, it is necessary to postulate the 
existence of new superfields carrying PQ charge. Although this extension of the MSSM 



1 Here G?„ is the gluon field strength tensor and Q a ^ v its dual. 
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can be realized in a number of ways, we will call the resulting weak-scale effective theory 
the PQMSSM. 

It has been noticed early on [11] that the PQ scale f a falls within the desired range 
for the SUSY breaking scale {msusy) in gravity-mediated SUSY breaking models: 

m 3/2 ~ m 2 SUSY /M P ~ 1 TeV => m S usY ~ 10 11 GeV ~ f a , (1.3) 

where is the gravitino mass and Mp is the reduced Planck mass. As a result, several 
models have been proposed [11-13] to connect the SUSY and PQ breaking scales. This can 
be achieved if the axion supermultiplet has tree level interactions with the hidden sector 
responsible for breaking SUSY. 

However, once a grand unified theory is assumed, the U(1)pq symmetry can appear 
as an accidental global symmetry of the theory, as naturally occurs in several SUSY 
GUTS [14]. In this case, f a will naturally be of order Mqut and will strongly violate 
its 10 12 GeV upper limit. It is possible to protect f a from obtaining Mqut contributions, 
either by breaking the PQ symmetry at a lower scale or artificially suppressing the vac- 
uum expectation value of the axion supermultiplet. Nonetheless such mechanisms always 
require the introduction of new superfields or fine-tuned parameters only for this purpose. 

String theory has emerged as an attractive ultraviolet complete theory which can easily 
incorporate axion-like fields as elements of anti-symmetric tensors [15]. Many of the would- 
be axions become massive, while the remaining light fields obey a global PQ symmetry. 
A survey of a variety of string models [16, 17] indicates that while PQ symmetry is easy 
to generate in string theory, the associated PQ scale tends to occur at or near the GUT 
scale rather than some much lower intermediate scale. This is in apparent conflict with the 
simple limits on f a from overproduction of dark matter as discussed above. 

One solution to the apparent conflict which allows for f a ~ Mqut is to invoke a 
tiny initial axion mis-alignment angle 0% ~ 0.003. In this case, one must accept a highly 
fine-tuned initial parameter which might emerge anthropically. 

An alternative solution was proposed in one of the original papers calculating the 
cosmic abundance of relic axions [18]: perhaps additional massive fields are present in the 
theory, whose late decays can inject substantial entropy into the universe at times after 
axion oscillations begin, but before BBN starts. In Ref. [18], it was proposed that the 
gravitino might play such a role. Several subsequent works have also explored the issue of 
dilution of (quasi)-stable relics via entropy injection [19-26]. 

Therefore it is of interest to investigate under which conditions the PQ scale can be 
extended to the GUT scale, while avoiding the known experimental constraints. While this 
possibility has been suggested before [21,23], we wish to explore this case using detailed 
calculations of particle production rates coupled with recent constraints arising from Big 
Bang nucleosynthesis (BBN). Here, we investigate the implications of f a ~ Mqut from a 
phenomenological point of view. In order to keep our conclusions as general as possible we 
will avoid choosing a specific GUT theory or PQMSSM model whenever possible. 

In Sec. |2[ we will present the general features of PQMSSM cosmology used in our 
analysis and review several well known results for this model. We will then discuss the 
f a ~ Mgut scenario in Sec. || and the case of a light axino LSP in Sec. ||. We will 
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present the main difficulties associated with large f a values and show how they can be 
avoided in the PQMSSM framework. In Sec. ||, we discuss the scenario where rria ~ ^3/27 
where the neutralino is assumed to be LSP. In several respects, this scenario is more 
appealing than the light axino LSP case. Section ^ summarizes our main results and 
discusses the implications of unifying the PQ and GUT scales. The appendices contain 
explicit formulae for evaluating the axion and neutralino relic abundances in radiation-, 
matter- and decaying-particle- dominated universes. 



2. PQMSSM Phenomenology 

In order to implement the PQ mechanism in supersymmetric theories, PQ charges have to 
be assigned to the MSSM fields and new PQ superfields must be introduced. The axion 
superfield is obtained from linear combinations of other elementary (non-MSSM) fields 
and is a singlet under the MSSM gauge group. Even though the full field content of the 
PQMSSM is highly model dependent, it must contain an axion supermultiplet composed of 
a complex scalar field (</>) and a Majorana fermion (a). The complex scalar field is usually 
divided into its axion (a) and saxion (s) components: 

s + ia 



" V2 

and the fermionic component is named axino. 

Since the axion field is the U(1)pq pseudo-Goldstone boson, it is massless, except for 
anomalous corrections coming from the QCD chiral anomaly. For temperatures well above 
A (the QCD chiral breaking scale), the axion is essentially massless, while for T <C A 
the QCD chiral anomaly induces a non-zero mass for the axion field. The temperature 
dependent axion mass is given by [27,28]: 



m° a , if T < A 

0.018 (^) 4 , if A < T ' 



where m® = 6.2 x 10~ 3 GeV// a , A = 200 MeV and T always refers to the thermal bath 
temperature. 

In order to solve the strong CP problem, the axion field must have effective couplings 
to the SU(3) gauge fields of the form shown in Eq. ([L . 1|) . Although other non- minimal 
interactions are possible, they are strongly model dependent and will be neglected here. 

The supersymmetric version of Eq. fll.ip implies the following couplings for the saxion 
and axino fields: 

Ceff = ^ S -^( G % GaflV + *&fDH0) + ^^rW, TlSGW + O(a^). (2.2) 
In most models, the axino also couples to the U(l)y gauge boson and gaugino: 

r a ayCa^aix) , . 
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where C a yy is a model dependent constant of order 1. Couplings between the saxion 
and axion, as well as between axinos and fermions-sfermions can also exist, but are model 
dependent. 

If supersymmetry is unbroken, both the axino and saxion are degenerate with the 
axion field, hence massless, except for the tiny QCD anomaly contribution. However, once 
SUSY is broken the saxion field will receive a soft mass of order 7713/2- On the other hand, 
being the fermion component of a chiral superfield, the axino remains massless at tree level. 
Nonetheless, the axino can receive loop corrections to its mass [29] of order m^ 2 / fa ~ 1 
eV, for f a ~ Mgut- In this case the axino is the lightest super symmetric particle (LSP) 
in the PQMSSM model. However, depending on the PQMSSM model, m a can also receive 
0(1713/2) contributions from the supergravity potential [30]. In this case, the axino might 
remain as the LSP, or the lightest neutralino could be the LSP [31,32], with Tn% < m a . 
The case of a gravitino as LSP is considered in Ref. [33]. 

In the case of axino LSP, the Z\ can be long-lived due to its l/f a suppressed cou- 
plings to the axino; this has important cosmological implications. Assuming the minimal 
interactions of Eq's. ( |2.2j ) and we have the following decay rates: 

(1)2 2 



F(Zi 5 + 7 ) = a Y C2 aYY COS 2 d «4 1)2 m 3 (1 _ ™Lf ( 2 . 4 ) 
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where vjp = (Zi\B) is the bino component of the neutralino field in the notation of 
Ref. [10]. If the neutralino is the LSP, the axino will be long-lived instead. The decay 
rates for a — > Z\ +^/Z are given by the above equations with m a . However, if the 

axino is heavier than neutralinos or gluinos, new decay modes a —¥ Z^, ZiZ or gg open 
up. Since the decay rates for these are discussed in Ref. [32], we do not reproduce them 
here. 



Due to the first interaction term in Eq. (2.2), the saxion decay width to gg is model 
independent, with decay width given by: 

2 3 

rv 777 

Saxions also decay to gluino pairs, but its width is always well below that to gluon pairs, 
as shown in Ref. [34]. The saxion might also decay directly into two axions. The decay 
width to axion pairs is given by: 

r(a -> aa) = i^m 3 s (2.6) 
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where k is a model dependent coupling. For our subsequent analysis, we will assume that 
the above decay mode is suppressed with respect to the one in Eq. ( |2.5[ ); this suppression 
is common in models with universal soft SUSY breaking terms [35]. 

2.1 PQMSSM cosmology 

The cosmology of PQMSSM models is very rich and has to be carefully examined. Here, 
we will assume that the PQ symmetry breaks before the end of inflation so as to avoid 
domain wall problems [36]. 

Due to their suppressed interactions, the axion, axino and saxion rapidly decouple 
from the thermal bath in the early universe or are already produced out of equilibrium, if 
the reheat temperature after inflation is smaller than the decoupling temperature [37]: 

10 12 GeV J V a 



'** = 10 GeV 77712^7 — ■ ( 2 - 7 ) 



However, for the values of f a considered here (/„ > 10 13 GeV), the decoupling temperature 
is always above f a . Since we only consider the case where the PQ symmetry is broken 
before inflation ends (Tr < f a ), we always have Tr < Td cp and axions, saxions and axinos 
are never in thermal equilibrium. In this case, the axion, saxion and axino thermal yields 
at T <C Tr are estimated as [35,38-40]: 

yTP 



g s J VlO 14 GeVy V fa 
T R \ /10 12 GeV\ 2 



10 14 GeVy V /, 



2 



where g s is the strong coupling constant at T = Tr and we have used <7*(T^ cp ) = 5*(Tr) = 
229. 

In an analogous way, gravitinos are thermally produced in the early universe with yield 
given by [41]: 



where & = (0.653, 1.604, 4.276) xlO^ 12 , h = (1.266,1.312,1.271), g t are the gauge couplings 
evaluated at Q = Tr and Mi are the gaugino masses also evaluated at Q = Tr. To compute 
the gravitino yield we assume mg > Mi, since we only consider cases with > 1 TeV. 
However, most of our results are weakly dependent on this assumption. If entropy is always 
conserved from Tr to To, Eq's. (|2.8j ) and (|2.9|) are still valid at T = To and can be used to 
compute the relic energy density today. 

Besides being produced from scattering of particles in the thermal bath, the saxion and 
axion fields can also contribute to the energy density through coherent oscillations. For 
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T S> m s (m a ), the expansion rate of the universe is too large and suppresses any oscillations, 
which only start at 3H ~ m s (m a ). The axion and saxion oscillation temperatures are then 
defined by [42,43]: 

3H(T a ) = m a (T a ) and 3H(T S ) = m s (2.10) 

where m a (T) is the temperature dependent axion mass. In most regions of the PQMSSM 
parameter space, we have T a ~ 1 GeV and T s ~ 10 10 GeV for a TeV scale saxion. 2 For 
reheat temperatures Tr smaller than T s , the saxion field starts to oscillate during the 
inflaton-dominated universe and will be diluted due to the inflaton entropy injection. In 
this case, the saxion coherent oscillation density increases with Tr until Tr > T s , where it 
then becomes Tr independent. The specific details of the transition from the static to the 
oscillating regime are strongly model dependent, but can be parametrized by an arbitrary 
initial field amplitude. If the axion oscillation starts in a radiation dominated universe, for 
T <C T 0jS , the axion and saxion coherent oscillation Yields are given by 3 : 

4.8 x 10- 9 /(fi)^ 6/12 (r«) (totA^v) 776 , if A < T a 

3.4 X KT 11 /^)^ 7 (T a ) ( 10 i/ n Gev ) > if T « < A 



7.8 x 10- 5 GeV (l^) (tot^cv) (f ) , if T. < T fi 
/•9X10- 8 GeV (jfe) (tot^v) ' (f) ' , if T R < T s 



(2.12) 



where [28] = 



In 



7/6 

and Sj and 9if a are the initial saxion and axion field 



amplitudes and are expected to be of order f a . 

From the above expressions we see that for most purposes the PQMSSM parameter 
space can be restricted to: 

{f a , m a , m s , Si, Oi, T R } + SUSY spectrum. (2.13) 

For simplicity, we take C a YY = 8/3 as in the DFSZ [9] axion model or in the KSVZ [8] 
model with PQ quark charges 2/3. 

2.1.1 Early saxion dominated universe 

Due to the small couplings of the axion supermultiplet, its components decouple from 
the thermal bath at very high temperatures. For temperatures below m s for thermally 
produced saxions or below T s for coherent oscillating saxions, the energy density of the 
saxion field scales as R~ 3 . Since the radiation energy density scales as R~ 4 , if the saxion s 
is sufficiently long-lived, it will dominate the energy density of the universe. Defining T e and 



2 For saxion oscillation in a radiation dominated universe, T s = . 1o1/2 ™s m p _ 

V t^I' 2 (T S ) 

3 Our saxion Yield expressions differ from Ref. [35] by a numerical factor because they assume the saxion 
field begins to oscillate at H — m s while we take 3-ff = m s to be consistent with the usual axion oscillation 
condition. 
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Tjj as the temperatures at which the saxion dominated era starts and ends, respectively, 
we have: 

Ps (T e ) = p 7 (T e ) T e = Un s Y s (2.14) 



T s = H(T D ) ^T D = - ^^ 1/4 (2.15) 



where T s is the saxion decay rate given by Eq. ( BJp and Mp is the reduced Planck mass. 
If To > T e , the saxion decays before dominating the energy density and its effects can be 
safely neglected 4 . However, if To < T e , the universe becomes matter (saxion) dominated 
for To < T < T e . At T = To most of the saxions have decayed and the radiation dominated 
era resumes. 

Fig. H shows T e as a function of the reheat temperature for m s = 20 TeV, S{ = f a and 
f a = 10 15 , 10 16 and 10 17 GeV. As we can see, T e increases with Tr until T s < Tr, when 
the saxion starts to oscillate after inflation. In this regime, the saxion Yield and hence T e 
becomes independent of Tp. From Fig. [l], we see that for /„ ~ 10 16 GeV and m s = 20 TeV, 
T e can be as large as 10 6 GeV. 

In Fig. ^, we plot Tp as a function of m s for f a = 10 15 , 10 16 and 10 17 GeV. We see 
that To can span a wide range of values, with To < 0.1 MeV for a sub-TeV saxion and 
fa ^ 10 15 GeV. However, if To < 5 MeV, entropy will be injected during the neutron 
freeze-out and the neutron-proton ratio will be significantly diluted. Since the success of 
BBN predictions strongly constrains the neutron-proton ratio to the value obtained in a 
radiation dominated universe, we must have To > 5 MeV in order to have a neutron freeze- 
out during the radiation dominated era, thus preserving the successful BBN predictions. 
As shown in Fig. g this requires m s > 50 TeV for /„ ~ 10 16 GeV. This bound can be 
softened by including additional model-dependent saxion decays such as s — > aa. 



4 If the saxion field has a large branching ratio into axions, its decays contribute to a hot DM component 
even for T e < Tr>- However, as mentioned in Sec. here we assume BR(s — > gg) ~ 1. 
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Once s starts to decay, entropy will be injected and effectively dilute all other relic 
densities already decoupled from the thermal bath, such as axinos, axions, gravitinos and 
possibly neutralinos. If the decoupling happens before the saxion dominated era (this 
is always the case for thermal axions, axinos and gravitinos), the dilution factor (r) is 
approximately given by [44]: 

Sf 1/4 Y q m q Tp 

T -i~^wkfa m £- <216) 

where g* is the number of relativistic degrees of freedom averaged over the saxion decay 
period, which we approximate by g*(Tu). Fig. || shows contours of r values in the Tr vs f a 
plane. We see that, for Sj = /„ and f a ~ 10 16 GeV, values of r larger than 10 are easily 
obtained. 

However, if the decoupling from the thermal bath happens during the saxion dominated 
period, as it might be the case for neutralinos or for the beginning of axion oscillation, the 
dilution factor will be modified due to the faster expansion rate during the saxion dominated 
phase. As a result, the neutralino freeze-out temperature (Tf r ) and the axion oscillation 
temperature (T a ) will be modified, if Tr> < T a j r < T e . 

Finally, if the decoupling (or axion oscillation) happens after the saxion decay, there 
will be no dilution. In summary, for the case of an early saxion dominated universe (T e > 
Tjj), we have: 

' n* D /r , if T e < T x 



X 

&x > ^ Tdd <T X <T e 
Q° D , if T D < T x < T DD 

, if t x < t d 



(2.17) 



where Q x represents the axion, axino, gravitino or neutralino relic density, Q X D is the 
corresponding relic density in a radiation dominated (RD) universe (such as when Trj > T e ), 
T x is the decoupling or axion oscillation temperature and Trj d marks the transition between 
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Figure 3: Entropy dilution factor (r) contours in the f a -Tn plane for m s = 5 and 50 TeV. The 
curves have r values 10 4 , 10 5 and 10 6 fr om bottom to top. The dashed region is excluded by the 
BBN constraints on Tr>. 



the matter dominated phase (MD) and the decaying particle dominated phase (DD) and 
can be approximated by (see Appendix): 

T DD = ( 9 -^T e T D . (2.18) 



g*(T x 

As mentioned before, the axino and gravitino fields will always have T x 3> T e , so 
their relic densities will simply be diluted by r. On the other hand the neutralino freeze- 
out temperature and the axion oscillation temperature may fall in any of the intervals in 



Eq. ( 2.17 ). The appropriate expressions for T a , Tf r , £l a and J)^ for the RD, MD and DD 



cases are given in the Appendix. To illustrate the behavior of and £l a in the distinct 
regions of Eq. (gl7b , we show in Fig. ^| the relic densities versus T e , with To = 5 MeV, 
9i = 0.4 and = 150 GeV, for f a = 10 15 , 10 16 and 10 17 GeV. As we can see, for small 
values of T e (< T a j r ), the axion and neutralino relic densities are simply diluted by 1/r. 
Once T e ~ 0.1 (10) GeV, the axion (neutralino) starts to oscillate (decouple) during the 
matter dominated (MD) era. As a result, the relic density is no longer diluted by 1/r, 
but by a smaller factor. Once T e 3> T a , the axion starts to oscillate during the decaying 
particle dominated phase (DD) and becomes independent of T e , despite the increase of 
entropy injection. Due to its large freeze-out temperature (Tf r ~ 7 GeV), the neutralino 
never decouples during the DD phase for the range of T e values shown. 



3. The f a = M GUT regime 

Assuming f a to be of order the GUT scale (~ 10 16 GeV) has several important consequences 
for the PQMSSM cosmology: 
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Figure 4: The axion and neutralino relic density versus T e , the temperature at which the universe 
becomes saxion dominated. The saxion decay temperature is fixed at its minimum value allowed by 
BBN (5 MeV), = 150 GeV and 9i — 0.4. The dotted, solid and dashed blue lines corresponds 
to the axion relic density for f a = 10 15 , 10 16 and 10 17 GeV, respectively. We also show the regions 
where the axion (neutralino) starts to oscillate (decouple) during the radiation (RD), matter (MD) 
or decaying particle (DD) dominated era. The neutralino relic density is normalized by its MSSM 
value. 



• The thermal production of all the components of the axion supermultiplet will be 
strongly suppressed by the large value of f a (see Eq. (2. 8)). 5 Therefore the axino 
and thermal axion contributions to the DM density are negligible and the saxion and 
axion densities are dominated by the coherent oscillation component. 

• Assuming ~ f a , large f a will result in a large energy density for the coherent 
oscillating saxion field as shown in Eq. ( |2,12| ), 

• Since T e oc f% (for coherent oscillating saxions) and Tp oc l/f a , large f a will usually 
result in a large entropy injection r = ^ from saxion decays, as shown in Figs. || 
and |. 

• The axion field will be extremely light (~ 10 -10 eV). 

• Since T s oc , saxions will be long-lived, potentially spoiling the BBN predictions. 

• For an axino (neutralino) LSP, ^z 1 (a) x V/a anc ^ the neutralino (axino) will be 
long-lived and a threat to successful BBN. 

From the above points, we see that f a ~ Mqut naturally leads to a long-lived saxion 
with large energy density from coherent saxion oscillations. As shown in Sec. 2.1.1, this 



5 The thermal production could still be relevant for Tr > f a . However, for such large Tr values the PQ 
symmetry would only be broken after inflation, resulting in a universe broken into domains with different 
values of 8i and Sj. Since this scenario requires a different dark matter treatment, we assume Tr < f a . 
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results in an early saxion dominated universe with a large dilution of other relics, due to the 
entropy injection during saxion decays. Below we discuss the implications of this scenario 
for big-bang nucleosythesis and the dark matter relic density. 

3.1 The dark matter constraint 

If the saxion field is neglected, the axion relic density is given by 



n a h 2 



9.23 x io-^ m) _l_ (.^L y , if T a < A 



\ 3/2 



1-32 02/(0,)-^, ( Tgri L) , if T a > A 



(3.1) 



where A = 200 MeV. Since the mis-alignment angle 9{ is supposed to be a random variable 
in the interval [0,7r], it is usually assumed that 8i ~ 1. Thus, the constraint on the dark 
matter relic density 

n DM h 2 = 0.1123 ± 0.0035 ( at 68% CL) (3.2) 

implies f a < 10 12 GeV. In this context the PQ and GUT scales can only be unified if 6i 
takes unnaturally small values (< 10~ 3 ). 

However, as shown in the Appendix, once the entropy injection from saxion decays is 
included, we have (for T a > A): 

' f a /Q x T D /T e , if T e < T a 
tt a <x I / a 14/n x T D /T e 4/11 , if T DD <T a <T e (3.3) 
j! /2 xT 2 , if T D <T a <T DD . 

Since Try oc l// a and T e oc / 2 (for coherent oscillating saxions), from Eq. ( |3.3[ ) we see that 
the axion relic density actually decreases with f a . In this case, the f a ^ 10 12 GeV bound 
can be potentially avoided. On the other hand, if T a < A, $7 a increases with f a , unless 
T a > T e . In both cases we see that £l a h 2 is maximally suppressed for Tp = T™ n = 5 MeV. 
From Fig. |], we see that the maximum dilution occurs in the DD regime with T a > A. Using 
the expression for £l® D h 2 in the Appendix, we estimate the maximum 0i value allowed by 
the DM constraint as: 

Hf»» < 0.11 * * < 0- = (69 - 106) (^f) ( 10 " /a GeV ) 3/4 . (3.4) 

The uncertainty on 8™' ax comes from the uncertainty on the axion mass at T a ~ A (see 
Eq. (O)). 



3.2 BBN bounds 

Since we will first assume a PQMSSM with a neutralino NLSP and an axino LSP, the 
neutralino will decay into axinos and SM particles and can be long-lived for f a ~ Mqut- 
Furthermore, the PQMSSM will also contain long-lived gravitinos and saxions. All or 
any of these three fields can decay during or after BBN, potentially spoiling its successful 
predictions, unless their relic densities are extremely small at the time of their decay. 
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The Z\ decay width and hadronic branching fraction are calculated in Ref. [45]. Fig. || 
shows the neutralino lifetime for a bino-like Zi, as a function of Tn% , for f a = 10 15 , 10 16 
and 10 17 GeV and mg Cm^. Usual BBN bounds on late decaying particles require [46] 

r < 0.01s if nh 2 > 1 or VLh 2 < 10 -4 if r > 10 4 s. (3.5) 

From Fig. |j| we see that unless is in the multi-TeV range, the neutralino life-time 
(t^ ) will be well above 10 4 s. To maintain sub-TeV values of m^, extremely small values 
of ^^/i 2 are required in order to satisfy the BBN constraints. Since in almost all of the 
MSSM parameter space 10~ 3 < 9, s i d h 2 < 10 3 [47], the BBN constraints would require an 
enormous fine-tuning of the MSSM parameters. However, if neutralinos decouple from the 
thermal bath before saxions have decayed, their relic density will also be diluted by the 



saxion decay, according to Eq. ( 2.17 ). As seen in Fig. [|, the neutralino dilution can exceed 
10 5 for large enough T e . Hence, the BBN bounds on late Z\ decays can be potentially 
avoided due to the large suppression of Q^/i 2 , without the need for fine-tuned MSSM 
parameters. 

Thermal gravitinos are produced out of equilibrium via radiation off of particles in the 
thermal bath (see Eq. (|2.9|) ) and have decay rates suppressed by 1/Mp, decaying during 



or after BBN, for mg < 50 TeV [48]. 6 Therefore, if Tr is large enough to significantly 
produce gravitinos in the early universe, their late decay will spoil the BBN predictions. 
This is usually known as the gravitino problem [49], common in most supergravity models 
with moderate and reheat temperatures above 10 5 GeV [50]. However, since gravitinos 
always have decoupling temperatures larger than the reheat temperature, the gravitino 
relic density is diluted by 1/r, as shown in Eq. ( 2.17f) . Since r can easily exceed 10 4 for 



f a ~ Mqut, as seen in Fig. y, the gravitino relic density will be strongly suppressed, 
naturally avoiding the gravitino problem. 



6 Here, we take mg to represent the physical gravitino mass, while m 3/ / 2 is the Lagrangian gravitino mass 
parameter. 
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In the scenario where m a ~ ""1.3/2 an d the neutralino is the LSP, axinos will cascade 
decay to neutralinos, as discussed in Ref. [32]. In this case the BBN bounds on late decaying 
axinos are easily avoided since Y? p is suppressed for large f a . Furthermore, if 



rria > m. 



91 



the decay mode a —> gg considerably reduces the axino lifetime and it usually decays before 
BBN. 

Finally, as already discussed in Sec. |2.1.1| , the BBN bounds on late decaying saxions 
require Td > 5 MeV. Using Eq's (2^5) and ( |2.15|) , we have: 



( f \ 2/3 

T D > 5 MeV => m s > 0.1 TeV ^ — ) . (3.6) 

Thus, as already shown by Fig. [2|, BBN bounds on saxion decays require m s > 50 TeV, for 
f a ~ 10 16 GeV. 

4. Results for an axino LSP 

In this Section, we discuss under which values of PQMSSM parameters we can conciliate 
f a ~ Mqut with the dark matter and BBN constraints for the case of a light axino as LSP 
and present explicit examples for our previous discussion. In this scenario, for large f a , the 
total dark matter relic density will be given by [51]: 



n DM h 2 = n a h 2 + nl p h 2 + ^% h 2 + ^n^h 2 (4.1) 



where we have neglected the subdominant contributions from thermal axions and the grav- 
itino contribution to the axion relic density, which is suppressed by m a /m^. To compute 
the above relic densities, we use Eqs. (PD , (|2l2|) , ( plTD and the expressions in the Ap- 



pendix. The BBN bounds on late decaying neutralinos and gravitinos are obtained from 
Figs. 9 and 10 in Ref. [46], with a linear interpolation for different values of or mg. 
For late decaying saxions we simply require To > 5 MeV, to avoid further entropy injection 
during the neutron freeze-out. 

4.1 Specific example for axino as LSP 

Fig. H shows the axion, axino, neutralino and gravitino relic densities as a function of f a - For 
the PQMSSM parameters we take m s = 50 TeV, s { = 10 f a , tuq = 1 TeV, T R = 10 11 GeV 
and m a = 0.1 MeV. We assume = 150 GeV and the neutralino relic density before 
dilution to be Q s ~ d h 2 = 10. For each f a value, a different value for 9{ is chosen so that 
JI^m/i 2 = 0.1123 is satisfied. As we can see from Fig. ^, the dilution of the axino, neutralino 
and gravitino relic densities rapidly increases with f a due to the increasing rate of saxion 
production via oscillations. For f a < 7 x 10 14 GeV, we have Slj?^ 2 ^ 10~ 4 and BBN 
constraints on late decaying neutralinos exclude this region. However, if a smaller value 
of £l s ~ d had been chosen, smaller f a values would be allowed. Once f a ~ 10 15 GeV, the 

Zi 

entropy injection from saxion decays dilutes the neutralino relic density to values below 
10 -4 , making these high f a values consistent with BBN. Finally, when f a ~ 10 16 GeV, 
the saxion starts to decay at T = Trj < 5 MeV and these solutions become once again 
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m~ = 150 GeV, £1- h 2 = 10, m, = l TeV, m = 50 TeV, m_ = 700 keV 




Figure 6: The axion, axino, neutralino and gravitino relic densities as a function of f a for — 
150 GeV, n^h 2 = 10, m a = 100 keV, m d = 1 TeV, m s = 50 TeV, T R = 10 11 GeV and s 2 = 10 f a . 
The misalignment angle (Oi) is chosen such as SloMh 2 ~ Sl a h 2 = 0.1123. The dashed region is 
excluded by BBN bounds on neutralino and gravitino late decays (for f a < 7 x 10 14 GeV) or 
To < 5 MeV (for f a > 10 16 GeV). 



excluded by the BBN constraints. We can also see that the gravitino relic density is 
strongly suppressed despite the large Tr value, easily avoiding the BBN constraints on 
late decaying gravitinos. Also, despite being the LSP, the axino does not significantly 
contribute to QoMh 2 , and the cosmologically allowed region around f a ~ 10 16 GeV has 
little dependence on ma- 
in Fig. we show the values of 9j necessary to satistfy the dark matter constraint for 
the same PQMSSM parameter values used in Fig. |6[ For f a < 10 14 GeV, the axion oscillates 
after the saxion has decayed (T a < Tjj) and fl a h 2 is not diluted by the early entropy 
injection. In this regime, the values of required to satisfy Q a h 2 = 0.1123 rapidly decrease 
with f a , since the axion relic density increases with f a for T a < Tjj. For 10 14 GeV < f a < 
6 x 10 15 GeV, the axion starts to oscillate in the decaying particle dominated (DD) regime 
(To < T a < Tdd) and il a h 2 decreases with f a , as discussed in Sec. 3.1. As a result, 8i 
increases with f a , although it is still required to be small (< 0.07). Once f a > 6 x 10 15 GeV, 
the axion oscillation still starts in the DD era, but now with T a < A. As shown in the 
Appendix, in this case Q a h 2 oc f a and the mis-alignment angle once again decreases as f a 
increases, although with a smaller slope than in the RD era. 

From Figs. || and f|, we see that f a ~ Mqijt can indeed be consistent with the dark mat- 
ter and BBN bounds. However, for the above choice of PQMSSM parameters, the region of 
parameter space consistent with all bounds is considerably restricted. Furthermore, B% still 
has to take small values, as would be the case in the PQ standard model cosmology, where 
the saxion field is neglected and f a ~ Mqut can be obtained if we take fl;<3x 10~ 3 [28]. 
Since the main purpose of the PQ mechanism is to avoid a huge fine-tuning in Oqcd, it is 
desirable to avoid unnaturally small values for the mis-alignment angle as well. With this 
is mind, we point out that 6i can take considerably larger values (~ 0.07) once the saxion 
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Figure 7: The misalignment angle required for ^luMh 2 = 0.1123 as a function of f a for — 
150 GeV, n^h 2 = 10, m a = 100 keV, m d = 1 TeV, m s = 50 TeV, T R = 10 11 GeV and s 2 = 10 f a . 
The dashed region is excluded by BBN bounds on neutralino and gravitino late decays (for f a < 
7 x 10 14 GeV) or T D < 5 MeV (for /„ > 10 16 GeV). 



and axino fields are included. Furthermore, the dilution of the neutralino and gravitino 
relic densities allows for an elegant way of avoiding the BBN constraints without having 
to assume extremely small Slg h 2 , low reheat temperatures or a multi-TeV gravitino. 

4.2 General results for axino LSP 

The above arguments are, however, limited by our choice of the PQMSSM parameters used 
in Figs. |6] and 0. In order to generalize these results, we perform a random scan over the 
following parameters: 

f a G [TO 15 , 10 17 ] GeV, 
m s e [10 3 , 10 5 ] GeV , 
Si/ f a G [10~ 2 , 10 2 ] , (4.2) 
Tr G [10 4 , f a ] GeV, 

and take = 1 TeV, m~ a = 100 keV, m ? =150 GeV and Q,% d h 2 = 10, as before. Our 
results will hardly depend on reasonable variation of these latter parameters, due to the 
enormous suppression of long-lived relics due to saxion production and decay. For each 
set of PQMSSM values the mis-alignment angle is chosen to enforce Q,DMh 2 = 0.1123. 
The BBN bounds on late decaying saxions, neutralinos and gravitinos are once again 
applied and solutions which satisfy all constraints are represented by blue dots. In order 
to differentiate the solutions excluded due to late decaying neutralinos or gravitinos from 
solutions excluded due to late decaying saxions (To < 5 MeV), we represent the former by 
red dots and the latter by green dots. 

Fig. [8| shows the scan results for the misalignment angle B% versus f a . As we can 
see, f a ~ 10 16 requires Oi < 0.07. Although small Oi values are still required in order to 
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Figure 8: The mis-alignment angle versus / a for a scan over the PQMSSM parameters as discussed 
in the text. All points satisfy J7dm^ 2 = 0.1123, while the green points are excluded by BBN 
constraints on late entropy injection from saxion decays (To < 5 MeV). The red points are excluded 
by BBN bounds on late decaying neutralinos or gravitinos and the blue points satisfy both the BBN 
and dark matter constraints. 



obtain f a ~ Mgut, the mis- alignment angle can now be twenty times larger than in the 
non-SUSY PQ scenario where the saxion/axino fields are neglected. Furthermore, if we 
require f a ~ 10 15 GeV instead, 6i can be as large as 0.4. We also point out that these 
conclusions are independent of our choice of ms and £l s ~ d , since the upper limit on Q, L 
comes entirely from the Td > 5 MeV constraint. These results also verify our estimate for 
gmax in Eq q ? wMch giveg g. < Q m _ .1(0.4 - 0.6) for f a = 10 16 (10 15 ) GeV. 

In Fig. ||, we show the saxion mass versus f a for the same points exhibited in Fig. |8|. 
As already discussed in Sec. |2.1.l| , the BBN constraint on late decaying saxions (Td > 
5 MeV) requires m s to be in the multi-TeV range, as shown by Eq. (|3.6]). Since we expect 
m s ~ ^3/2) models with large m 3 / 2 ~ 10 — 50 TeV such as Yukawa- unified SUSY [52], 
mirage unification [53,54], effective SUSY [55], AMSB [56] or string-motivated models such 
as G2-MSSM [57] would naturally yield such heavy saxions. 

We also see that the BBN bounds on late decaying neutralinos do not significantly 
constrain the saxion mass, since allowed (blue) solutions can be found for any m s value 
as long as the bound in Eq. ( |3.6[) is satisfied. We also point out that the few red points 
with saxion masses below the limit in Eq. (|3.6| ) correspond to the case T e < Td, where the 
saxions decay before dominating the energy density of the universe. All these solutions 
have extremely small 0j values, which lie in the narrow band at 9i ~ 0.001 seen in Fig. ||. 

Fig. |K] shows the saxion field amplitude s« versus f a . As discussed in Sec. Q s» 
parametrizes the details of the transition from the static to the oscillatory regime of the 
saxion field near T = T s . To compute the value of Sj, the full saxion potential for T >T S 
needs to be known, which requires assuming a specific PQMSSM model as well as knowledge 
of the SUSY breaking mechanism. Nonetheless, natural values for Sj are f a or Mp. Fig. 10 
shows that small values of Si/f a are disfavored, since they suppress the entropy dilution of 
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Figure 9: The saxion mass versus f a for a scan over the PQMSSM parameters as discussed in the 
text. All points satisfy VLoMh 2 — 0.1123, while the green points are excluded by BBN constraints 
on late entropy injection from saxion decays (To < 5 MeV). The red points are excluded by BBN 
bounds on late decaying neutralinos or gravitinos and the blue points satisfy both the BBN and 
dark matter constraints. 
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Figure 10: The initial saxion field amplitude (s^ divided by f a versus f a for a scan over the 
PQMSSM parameters as discussed in the text. All points satisfy ftoAih 2 — 0.1123, while the green 
points are excluded by BBN constraints on late entropy injection from saxion decays (To < 5 MeV). 
The red points are excluded by BBN bounds on late decaying neutralinos or gravitinos and the blue 
points satisfy both the BBN and dark matter constraints. 



the neutralino and gravitino relic densities, conflicting with the BBN bounds. However, as 

seen in Fig. 1C, unnaturally large Sj// a values are not necessary for obtaining f a ~ Mqut- 

Furthermore, smaller Q~ d h 2 values would allow for smaller Sj/ f a . 
i — i Zl 

Finally- in Fig. |11| - we show the reheat temperature versus f a . As already mentioned 
in Sec. |2.1.1| , large Tr increases T e , resulting in an increase in the dilution of the axion, 
neutralino and gravitino fields. From Fig. 11, we see that Tr > 10 8 GeV is usually required 
to satisfy the BBN constraints on late decaying neutralinos and gravitinos. We point out 
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that large reheat temperatures are motivated by thermal leptogenesis models [34, 45, 58] 
which usually require Tr > T R nm ~ 2 x 10 9 GeV in order to explain the observed matter- 
antimatter asymmetry (r) = nB n ns ~ 6 x 10~ 10 ) in the universe 7 . However, due to the 
entropy injection during saxion decays, the asymmetry hr — n B will also be diluted by 1/r. 
In this case the actual lower limit on the reheat temperature is Tr/t > T™ n . Hence, in 
frame b we show instead f a vs. Tr/t, and find that all BBN-allowed points have Tr/t < 
10 6 GeV. This rigid limit comes from the BBN bounds on late decaying gravitinos. 
For mg > Mi (as assumed here), where Mj are the gaugino masses, we have [48]: 

W - * x 10- (^) (£) and r s , tf. (1^)' ,4.3) 

where Q^h 2 is the gravitino relic abundance including the entropy dilution. Thus, for a 
1 TeV gravitino, we have rg ~ 10 5 s and the BBN bounds described in Sec. 3.2 require 



Og/i 2 < 10 -4 , which implies Tr/t < 10 6 GeV, as seen in Fig. lib. Therefore, the unification 
of the PQ and GUT scales seems to strongly disfavor thermal leptogenesis [59] scenarios, 
unless a heavier gravitino is assumed, as in the usual MSSM scenario. In particular, for 



= 10 TeV, Eq. (4.3) gives tp; ~ 10 2 s and the BBN bounds are considerably weaker 



G 

in this case [46]: Q^h 2 < 0.5. Hence, for multi-TeV gravitinos, we can have Tr/t < 
10 9 GeV, which makes thermal leptogenesis once again viable. We also note that non- 
thermal leptogenesis requires only Tr/v > 10 6 GeV [60], while Affleck-Dine leptogenesis 
allows still lower Tr/t values [61]. 

5. The neutralino LSP case: mixed aZ\ DM 

So far, our results have focussed on the case of the PQMSSM with an axino as LSP, so that 
dark matter consists of an ad mixture. Our results were largely independent of reasonable 
variations in rria since the axino abundance suffers a huge suppression due both to the 
large values of f a which are required and to entropy production from saxion decays. A 
qualitative difference results if we take m„ so high that > nt^ and the neutralino Z\ 
becomes the LSP, so dark matter would consist of an aZ\ mixture. In this case, gravitinos 
and axinos can still be produced thermally at high Tr, but now these states will cascade 
decay down to the stable Z\ state, and possibly add to the thermal neutralino abundance. 

In the mixed aZ\ DM scenario, neutralinos are produced via axino decay at temper- 
ature T D = y/TiMp I (ttV (T d ) / 90) V 4 , as well as via thermal freeze-out at T = Tf r . The 
neutralinos from axino decay may re- annihilate at T = Td if the annihilation rate exceeds 
the expansion rate at To [31]: 

s s 

where (Td) is the total neutralino number density due to thermal (freeze-out) and non- 
thermal (axino decays) production. Thus, the re-annihilation effect depends on a large 



7 Here, ub is the baryon density, rig is the antibaryon density and n 1 is the photon density of the 
universe. 
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Figure 11: In frame a, we plot the reheat temperature versus f a for a scan over the PQMSSM 
parameters as discussed in the text. All points satisfy floAih 2 = 0.1123, while the green points 
are excluded by BBN constraints on late entropy injection from saxion decays (Td < 5 MeV). The 
red points are excluded by BBN bounds on late decaying neutralinos or gravitinos and the blue 
points satisfy both the BBN and dark matter constraints. In frame b, we show the value of Tr/t 
to account for possible dilution of the baryon asymmetry due to entropy production from saxions. 



thermal production rate for axinos. In the f a ~ Mqut case considered here, neutralino 
re- annihilation is largely irrelevant because 1. thermal production of axinos is suppressed 
by 1/ f% and 2. the axino abundance at T = Tp is also highly suppressed by the saxion 
entropy production. Thus, for the large f a scenario, the neutralino abundance is estimated 
to be: 

% = n T / + + (5.2) 

1 z i rria m c 

where is evaluated for either a MD, DD or RD universe, and Hj p and f2~ p are diluted 
by entropy production ratio r for r > 1. At the end, we must add in the axion abundance 
f2 a as calculated for a MD, DD or RD universe, with the latter case diluted as usual by 
entropy ratio r when T e <T a . 

Our results in Sec. || showed that models with a very light axino as LSP could be 
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consistent with f a ~ Mqut, but only with very large- perhaps uncomfortably large- 
values of the saxion mass, with m s typically in the tens of TeV range. In gravity-mediated 
SUSY breaking models, a puzzle would then arise as to why the sparticles exist in the 
sub- TeV range, while saxions are present at 10-50 TeV. 

Here, we note that there do exist several SUSY models where ra s is naturally at the 
tens of TeV scale [52,55,56]. One possibility consists of models with mixed moduli-anomaly 
mediated SUSY breaking soft terms (mirage unification, or MU) [53]. These models are 
based on the KKLT proposal [62] of string models where the moduli fields are stabilized via 
fluxes and the addition of a non-super symmetric anti-D-brane breaks SUSY and provides 
an uplifting to the scalar potential leading to a deSitter vacuum. In this class of models, 
Choi and Jeong have calculated the magnitude of soft SUSY breaking terms where the 
strong CP problem is solved by the PQ mechanism [54]. Since the MSSM soft terms arise 
from mixed moduli/anomaly mediation, their magnitude is at the TeV scale even though 
7TI3/2 is naturally in the multi-TeV regime. Furthermore, they find that typically m„ ~ mg, 
and the saxion mass m s ~ \/2m^. 

Another possibility arises from M-theory models with compactification on manifolds 
of G2 holonomy. In these models [57], the gravitino, axino and saxion masses, along with 
MSSM scalars, are all expected to inhabit the multi-TeV range, while gaugino masses are 
expected to be much lighter, and non-universal. These latter models typically yield a 
wino-like lightest neutralino. 

The MU soft terms have been programmed into Isajet/Isasugra [63], and are functions 
of the mixed moduli-AMSB mixing parameter a, 7713/2, tan/3 and sign{ij). They also 
depend on the matter and Higgs field modular weights rtj, which can take values of 0, 1/2 
or 1, depending on if the fields live on a D3-brane, a D7 brane or an intersection. Spectra 
for the string models based on G2 holonomy can readily be generated with Isasugra using 
non-universal gaugino mass entries. 

In Fig. [l^a)., we show the axion and neutralino relic abundances versus f a for the 
mirage- unification model with moduli/AMSB mixing parameter a = 3, m 3/ / 2 = 40 TeV, 
tan/3 = 10 and \i > with m t = 173.3 GeV. We take m„ = mg = 40 TeV, m s = v2mg = 
56.6 TeV and Tr = 10 11 GeV. We further take modular weights n m for matter fields 
equal to 1/2 and for Higgs fields nn = 1- The Isajet spectra gives m g = 877.4 GeV and 
niy = 352 GeV and Q s ~ d h 2 = 0.026, where Z\ is mainly higgsino-like. From the plot, we 
see that for f a < 1.5 x 10 14 GeV, too much neutralino dark matter is produced due to 
axino decays and already saturates the DM constraint, even for f2 a /i 2 = 0. As f a increases, 
the axino abundance falls sharply since the thermal production rate is suppressed by 
and the Yield is diluted by entropy injection from saxion decays. The gravitino abundance 
also falls, but not as sharply, since here the diminution is only due to entropy dilution 
from saxion decays. The thermal neutralino abundance falls, but less sharply still, since, 
for the choice of parameters in Fig. |l^, the neutralino freezes-out in a MD universe, which 
results in a dilution smaller than 1/r, as shown in Fig. The relic axion abundance grows 
with f a , but is also diluted. Here, we dial 6i to an appropriate value such that h 2 is 
fixed at the mesured value of ~ 0.11. For f a > 4 x 10 GeV, the dark matter is axion- 
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dominated. Once we reach values of f a > 10 16 GeV, then the saxion decay temperature 
Td drops below 5 MeV and we consider the model BBN excluded (dashed curves). When 
compared to Fig. |6[ we see that the neutralino LSP case easily avoids the BBN bounds on 
late decaying axinos, as already mention in Sec. 3.2. 

In frame b.), we show the value of 9i which is needed for the case of m 3 / 2 = 40 and also 
for 30 TeV. The value of 9{ is again typically in the 0.04 — 0.08 range in order to suppress 
overproduction of axions. 

The main results from the scan over parameter space performed for the axino LSP 



case, shown in Fig's. B-ll, still hold for the neutralino LSP scenario, since they are weakly 



dependent on the nature of the LSP. However, since the BBN bounds on axino decays 
are easily avoided in the heavy axino case, the lower bounds on Tr and Sj// a are now 
relaxed in the neutralino LSP scenario. Nonetheless, the upper bound on Tr/V, relevant 
for baryogenesis mechanisms, still holds (if we keep mg = 1 TeV), since it only depends 
on the gravitino mass, as discussed in the last Section. 



6. Conclusions 

We have investigated the possibility of extending the usual upper bound on the Peccei- 
Quinn scale (f a ) to values of the order of the Grand Unification scale (Mqut)- We show 
that large f a values usually lead to an early universe dominated by coherent oscillating 
saxions, which are required to decay before Big Bang nucleosynthesis. In the case of a 
light axino as LSP, the injection of entropy during the decay of the saxion field results in 
a dilution of the axion relic density, which allows us to evade the usual upper bound on 
fa (^ 10 12 GeV). Furthermore, the dilution of the neutralino and gravitino relic densities 
naturally evade the BBN bounds on late decaying Z±s and Gs. From Eq's. ( |3.4j ) and (|3~G|), 
verified by the scan over parameter space, we find that, in order to allow f a ~ 10 16 GeV: 

• 6i < 0.07 — 0.1 is necessary to satisfy the axionic dark matter relic density constraint, 

• m s ^ 50 TeV in order to satisfy BBN constraints on late decaying saxions. 

Furthermore, for the axino LSP case with a neutralino LSP with = 150 GeV and 
h 2 = 10: 

• si/ f a > 1 and Tr > 10 s GeV are required to increase coherent saxion production 
and hence increase dilution of the neutralino relic densities and to satisfy the BBN 
bounds. 

While the first two conditions are quite independent of the SUSY spectrum chosen (parametrized 
here by tt s ~ d h? and m^), the third condition can be relaxed if PQMSSM models with heav- 
ier neutralinos and/or smaller Q s ~ d are considered. 

Z\ 

We also investigated the case where m a ~ 7713/2 an d the neutralino is LSP such that 
dark matter is comprised of an axion/neutralino mixture. Models such as Mirage Unifica- 
tion or string models based on G2 holonomy naturally give axino and saxion masses in the 
tens of TeV range, while maintaining at least some superpartners below the TeV scale. In 
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MM1: m- = m~ = 30 TeV, m s =\f2m- 




Figure 12: The axion, axino, neutralino and gravitino relic densities as a function of f a for the 
Mirage Unification model of SUSY breaking, with a — 3, m 3 / 2 = 40 TeV, tan (3 — 10 and [i > 0. We 
also take ma = 1713/2 and m s = \/2m^/2: with Tr = 10 11 GeV and s, = 10/ a . The misalignment 
angle (6i) is chosen such as f^j^ft 2 — Q a h 2 — 0.1123. The dashed region is excluded by BBN 
bounds on saxion decay: Tjj < 5 MeV. In frame b)., we show the required value of misalignment 
angle 9i for the case of 17*3/2 = 40 TeV and also for 777,3/2 = 30 TeV. 

these models, if f a ~ Mgut, then again we expect large amounts of entropy production 
from saxion decay, while neutralino, axino and gravitino abundances are all supressed to 
tiny levels, thus helping to avoid BBN constraints. 

While all our results were caluclated assuming the saxion s — > gg decay mode at 100%, 
we note that other model-dependent decay modes such as s — >• hh or s — > aa may be 
present. The first of these would contribute to additional entropy production and decrease 
the saxion lifetime, thus helping to avoid BBN constraints. In this sense, we regard our 
results as conservative. On the other hand, if saxion decay into axions is significant, then 
saxion decays inject relativistic axions and increase the effective value of g* during BBN. 
In addition, the entropy injection from saxion production and decay is greatly diminished. 
In this case, in order to keep the sucessful BBN predictions, p s has to be suppressed [35], 
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which likely disfavors the f a ~ Mqut scenario. 

As consequences of the f a ~ Mqut scenario, we would expect the DM of the universe 
to be axion-dominated, with a tiny component of either axinos or neutralinos. The axion 
mass is expected to lie in the 10 -10 eV range which is well below the range currently being 
explored by the ADMX experiment [64]. New ideas and new experiments will likely be 
needed to explore the axion direct detection signal in this mass range. Furthermore, the 
f a ~ Mqut case can accommodate a much wider range of ^1%^? values than the pure 
neutralino DM scenario in MSSM models, since, for a neutralino LSP, its relic abundance 
is suppressed by the entropy injection, while, for an axino LSP, the neutralino contribution 
to the DM relic abundance is suppressed both by entropy dilution and ms/m^r . 

Note added: As this manuscript was nearing completion, a paper by Kawasaki et al. 
Ref. [65] was released; they also show that large f a ~ Mqut can occur in an inflationary 
model which relates inflation to the PQ scale. 
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A. Relic Densities in a Saxion Dominated Universe 



Here we briefly review the cosmology of an early saxion dominated universe and present 
the expressions for the axion and neutralino relic densities used in Sec. |3[ 



As discussed in Sec. 2.1.1, we assume that the saxion field becomes the main component 



of the universe's energy density at a temperature T e given by Eq. (2.14). Therefore, the 
universe becomes matter dominated until the saxion decays at T = Tp>, with the decay 



temperature given by Eq. (|2.15| ). 

During the saxion dominated phase (To < T < T e ), the universe expands at a faster 
rate, given by: 

3/2 

(A.l) 



T 



(R(Te)\ 



M P V R(T) J 

where g*(T) counts the number of relativistic degrees of freedom 8 , R(T) is the scale factor, 
Mp is the reduced Planck Mass and we have used Eq. ( [2. 14 ). For T S> Tp>, most of the 
saxions have not yet decayed and entropy is conserved: 



g*(T e )Ti 



(R(Te)\ 

V R(T) J 



Therefore, 



7T 2 T 3 / 2 
H(T) = - x l— g JT)T e 



{T D <^T <T e ). 



(A.2) 



(A.3) 



For simplicity, we assume (?* (T) = g*s(T), where p rad = tt 2 g,(T)T 4 /30 is the radiation energy density 
and s = 2ty 2 g*s(T)T 3 /45 is the entropy density of the universe. 
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On the other hand, if T > Tp, the entropy injection of saxions results in [19]: 

<?*(T) 2 T 8 (R{T D 



9*(T D ) 2 T 8 D " 
H{T) = H{T D 



R(T) 
3* COT 4 



g*(T D )T*' 

Since at T = To the energy density is once again radiation dominated, we have: 



H(T D ) 
» H(T) 



7T 



^ 9*{T) 



(T D < T < T e 



(A.4) 



A.l The Axion Relic Density 

Due to the distinct expressions for H(T) during the matter dominated (MD) or decaying 
particle dominated (DD) phases, the axion oscillation temperature (T a ) as well as the co- 
herent axion oscillation relic density can significantly deviate from the usual expressions, if 
the axion starts to oscillate in the saxion dominated era (Td < T a < T e ). Using Eqs. ( |A.3| ), 
( |A.4j ) and Eq. ( |2.1| ) for the T dependent axion mass and Eq. ( 2.10 ) to define T a we find: 



a, 



f fiAD/r , if T e < T a 
n™ D , if T DD <T a <T e 
n° D , if T D < T a < T DD 

[n* D ,i£T a <T D 



(A.5) 



where r is defined in Eq. ( 2.16| ) and 9 

9.23 x 10- 3 ^ 2 /(^ 



n? D h 2 



fa 



1 / fa 



i.32 elm) 

7.5 x 10- 5 e 2 f[H)T L , 

i-4 efm) 



fa 

To 15 



fa 



7/6 
2 

14/11 



3/2 



, if T a < A 



if T a > A 



if T a < A 



ifa, if r a > A 



7.5 x 10- & 6if{6,i)T D v l0 



fa 



1.72 



rp2 [ fa 



, if T a < A 

3/2 



if T a > A 



(A.6) 



(A.7) 



(Ai 



The corresponding expressions for the oscillation temperature are: 

1/2 



1.23 x 10 2 



10 1 



8.71 x 10 



-l 



?.(T )Vi 



10 12 

/a 



,if^ D <A 
1/6 , if 1*° > A 



(A.9) 



9 All quantities are in GeV units. 
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iDD 



0.11 x 10 2 



9.0 x 1Q- 1 



2 ( V9*(Td) 



g,{T a ) f a 1 D 



1/4 



if T? D < A 



if T DD > A 



(A.ll) 



The temperature Tdd marks the transition from the matter dominated phase to the decay- 
ing particle dominated phase, where entropy is no longer conserved. An approximate value 
for Tdd can be obtained matching the axion relic densities in the MD and DD phases: 



Tdd 



g*(Tj 



D. 



9*(T a ) 



T e T D 



1/5 



(A.12) 



A. 2 The Neutralino Relic Density 

The neutralino will decouple from the thermal bath when 

(av)n Si (T fr ) = H(Tf r ), 



(A.13) 



where Tf r is the freeze-out temperature, (av) is the neutralino annihilation cross-section 
and 

(m~ T\ 
-fr) e ~ m ~ zJT - 

In a radiation dominated universe, the neutralino yield is given by: 

H(T fr ) 



Y z (T fr 

while in a matter dominated universe 

Y z (T fr ] 



(av)s(T fr y 

3 H(T fr ) 
2{av)s(T fr y 



(A.14) 



(A.15) 



(A.16) 



As in the axion case, the neutralino can freeze-out before the universe becomes matter 
dominated (Tf r > T e ), during the matter dominated phase (To <C Tf r < T e ), during 
the decay dominated phase (To < Tf r <C T e ) or during the radiation dominated phase 
(T fr <T D ). Using Eqs. flO}, flO} and ( |Al3| )-( [A~T6| ), we obtain: 



{ n^/r , if T e < T fr 
Vf D , if T DD < T fr < T e 
, if T D < T^ < T DD 

, if Tfr < T D 



(A.17) 



where 



Q^h 2 = 8.5 x IP' 11 - 1 ( ^ 

^ 2 yr^: * 



* 2Y g*(T fr )Tj r * 

The corresponding freeze-out temperatures are given by: 



/ GeV 2 

V M 



(A.18) 
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(A.19) 



Once again T^d marks the transition between the MD and DD phases and is given by: 



Tdd 



D) 



9*{ T fr) 



1/5 



(A.20) 
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